Anomalous Scattering of Low-lying Excitations in a Spin-1 Bose— Einstein Condensate 
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We present the simplest theory of perfect tunneling of an excitation in a Bose-Einstein condensate 
(BEC) through an impurity potential with an arbitrary shape in the low-momentum limit. That 
is for the transverse spin wave in the ferromagnetic phase of a spin-1 BEC. This mode obeys a 
Schrodinger-type equation; yet, effects of the potential on its transmission coefficient T and on 
its scattering cross section a vanish in that limit. The order parameter determines T, and the 
momentum p-dependence of a exhibits a Rayleigh scattering type (a tx p 4 ). These properties are 
common between two types of Nambu-Goldstone modes: this spin wave and the Bogoliubov mode. 
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Introduction. — Tunneling through a barrier without 
reflection is a curious phenomenon. Such anomalous 
behavior is observed in superfluids and superconduc- 
tors. In the early 2000s, perfect transmission in the 
low-momentum limit against a potential barrier was pre- 
dicted in [l], Q for the low-lying excitation (the so-called 
Bogoliubov mode Q) in a Bose-Einstein condensate 
(BEC). This perfect transmission was termed anomalous 
tunneling Q. Subsequently, many interesting tunneling 
properties have been discovered in scalar BECs [i - 11 1 . 

In quantum electrodynamics, it is known that a rela- 
tivistic particle tunnels through a high and wide barrier, 
which is referred as the Klein paradox (l^ |: it is also a 
recent subject of graphene [13[. This counterintuitive 
phenomenon is described by the Dirac equation for two- 
component wave functions (WFs) and its particle has the 
linear dispersion. Since the Bogoliubov mode has the lin- 
ear dispersion in the low-momentum limit and obeys the 
Bogoliubov equation for the two-component WFs, one 
might expect that this paradox is linked to anomalous 
tunneling of the Bogoliubov mode via some ways like the 
charge-conjugation symmetry for the Dirac equation. 

Two years ago, however, anomalous tunneling through 
a generic barrier was derived using the coincidence be- 
tween the WF of the Bogoliubov mode in the low-energy 
limit and that of the condensate [Ej], which is regarded 
as a property of NG modes. With this explanation, we 
find ourselves confronted with the following questions 
which we take up below: Is anomalous tunneling inherent 
to Nambu-Goldstone (NG) modes in broken-symmetry 
states? Does a NG mode with a parabolic dispersion 
relation described by a Schrodinger-type equation (cor- 
responding e.g. to the ferromagnetic magnon) exhibit 
anomalous scattering properties? 

The ferromagnetic state of the spin-1 spinor BEC [H], 
[l5l | (a BEC composed of particles with spin-1 internal 
degrees of freedom) is suitable for investigating these un- 
resolved questions since both gauge and spin rotational 
symmetries are spontaneously broken and accordingly 
the system hosts two types of NG modes: the Bogoliubov 
mode (type I) and the spin wave (type II) [3] • By com- 



paring the scattering properties of these two modes, we 
can gain a deeper understanding of the scattering prop- 
erties inherent to NG modes in broken-symmetry states. 
Further, the study of the spin wave in the spinor BEC 
is helpful for understanding the scattering properties of 
spin waves in magnetic systems; tunneling of spin waves 
in a magnetic film through an inhomogeneous magnetic 
field was recently demonstrated 17|, |l8[ . 

With this background, we investigate the scattering 
properties of low-lying excitations in the ferromagnetic 
phase of a spin-1 BEC from an obstacle within mean- 
field theory. We examine their tunneling properties and 
their cross section a of scattering from a generic potential 
with spherical symmetry. Our main findings on the spin 
waves are that anomalous tunneling occurs at the low- 
energy limit, and that the Rayleigh scattering type k- 
dependence a cx k A holds at small wavenumbers k. Com- 
parison with the Bogoliubov mode and inspection of the 
derivation suggest that both findings hold for other NG 
modes regardless of the spectrum. 

Formulation. — We start with a Hamiltonian for spin-1 
bosons with the atomic mass to, given by [lij 
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where \&j is the field annihilation operator for a spin-1 bo- 
son in hyperfine state j (j — ±1, 0), S a {a — x, y, z) are 
spin matrices, and repeated indices are summed. When 
we denote by cif the s-wave scattering length for the to- 
tal spin F channel, two coupling constants are given by 
Co = 47rft 2 (ao + 2d2)/(3m) and C\ = 4irh 2 (a2 — ao)/(3rn). 
The external potential V ex t is coupled only to the local 
density, so that H is U(l) x SO(3) invariant (l5| . 

The mean- field theory of the spin-1 BEC was estab- 
lished in 

0, 

[lij . The ground state of the ferromag- 
netic phase is realized for c\ < 0, and its condensate 
WFs (order parameters) (Wj) = <&je~ l f lt / h are given 
by ($ + i,$ ,$-i) = (vW>(r),0,0) [H, where fx is the 
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chemical potential and n is the density of the conden- 
sate in the uniform regime where </>(r) becomes unity. 
We note that <f> obeys an equation which is identical in 
form with the GP equation of a scalar BEC: Hgp4> = 0, 
with H G p = -{h 2 V 2 )/(2m) + V cxt + c + n(cj) 2 - 1) and 
c + = Co + ci, where we used fi = c + n [14]. In the ground 
state, we can take 4> as real. This system includes three 
types of excitations 3 HI- The excitation of j = +1 
is mathematically identical to the Bogoliubov mode of 
the scalar BEC; its tunneling properties can be drawn 
from results that have been established in previous stud- 
ies P, 0, 0" U|- The spectrum of the mode j = — 1 is 
massive: E = fi 2 P/(2m) + 2|ci|n. Since our emphasis is 
on the nature of NG modes, we will not discuss further 
this massive excitation in this Letter. Its properties have 
been partially reported in and will be reported else- 
where. The mode of j = leads to the transverse spin 
wave with the spectrum E = h 2 k 2 / (2m), which we focus 
on in this study. When cf>o denotes the WF of the exci- 
tation for j = 0, its equation of motion is given by [H[ 

ih— — ='Hgp0o, with ct> = (()oexp(—iEt/h). (2) 
ot 

We introduce dimensionless units: E = E/(c+n), V ext = 
Vext I (c+n) , k = fc£f, and x = x/£f where £f = hyfmc+n. 
We hereafter omit the bar for simplicity. 

Transmission coefficient. — We first study reflection 
and refraction of the transverse spin wave. We consider 
a potential Vext( r ) = V ex t(x) with only an x dependence, 
where V cx t(x) acts as a potential wall near x ~ 0. We 
further impose the following condition: V cx t(x) — s- Vr 
for x >■ 1, and V ext (x) — > Vl for x <C — 1. We as- 
sume that the incident wave with momentum k propa- 
gates at an angle #l with respect to the x axis [201 ] . 
For |af| — > oo, we find <j> — > \/l — V v (= A v ) through the 
equation of <j>, where A v is the amplitude of the order 
parameter |$ + i|(= |(^ + i)|), and the subscript v denotes 
L (or R) for the asymptotic regime of x < (or > 0). In 
dimensionless form, the healing length is modified such 
that CL(orR) = l/4(orR) for the regime x < (or > 0), 
since the condensate density is not unity at \x\ — > oo. 
Even if Vl ^= Vr, there is a common energy spectrum 
E = h 2 k 2 /(2m) for \x\ > £l,r from ©. From this result 
and the continuity of the WF, it follows that the spin 
wave is not refractive. This is to be contrasted with the 
behavior of the Bogoliubov mode, which was shown to 
be refractive in Q , obeying Snell's law. 

We take the xy plane as the incident plane without 
loss of generality, and set </>o( r ) = X(x)exp(iky sin 9^). 
Let k<£ — k cos #l be the x-component of the momentum, 
and let us define Hqp = + V ex t + <t> 2 - 1- For 

elucidating tunneling properties, we solve 



with the following boundary condition: 



H GP X(x) 



k 2 

fX(x), 



(3) 



X(x) 



exp (ik x x) + rexp (—ik x x) for x <C — £l 



t exp (ik x x) 



for x » £r. 



(4) 

We denoted by r = r(k x ) and t — t(k x ) respectively 
the amplitude reflection and transmission coefficients. 
The transmission and reflection coefficients are given by 
T(k) = \t(k)\ 2 and R(k) = \r{k)\ 2 . 

Below, we prove that these coefficients are given by 



lim T(k) 



4M L M B 



k^0 



(M L + M R ) 2 ' fe^o 



lim R(k) 



(M L - M R ) 
(Ml + Mr) 



(5) 



Here, M v for v — L and R denotes the modulus of 
the magnetization |M| in each asymptotic regime; it is 
expressed as [i^nA 2 when we recover the conventional 
units. This result satisfies T + R= 1, as well as the usual 
reciprocity relationship for tunneling. When Vl = Vr, 
the magnetizations on the left and right sides are equal 
(i.e., Ml = Mr) and perfect transmission at low energy 
T(k — > 0) = 1 follows regardless of the forms and heights 
of the barriers about x = 0. This is the simplest exam- 
ple of anomalous tunneling; earlier anomalous tunneling 
is derived through a two-component Schrodinger equa- 
tion for the Bogoliubov mode [l|, 0], while anomalous 
tunneling of this spin wave is derived through a single- 
component Schrodinger equation. We find that this tun- 
neling phenomenon is nontrivial by recalling that a par- 
ticle obeying the single- component Schrodinger equation 
generally undergoes total reflection against a potential 
barrier at k 0. 

Let us prove JS]). We first seek the solution to ^ in 
the form of X(x) = J2n=o kl n X^(x). X {n \x) satisfy 

H GP X^(x)=0, H GP X^(x) = ~X( n - 1 \x) (6) 

for n > 1. As a set of linearly independent solu- 
tions to the first equation of ([6]), we take X\(x) = 
4>(x) and X\\(x) = —j(f>(x) + <fr(x) dx' 4>~ 2 {x'). Here, 

7 = ( 7 R + 7l)/2, 7 l = f-°°dx(4,- 2 (x)-A^), and 
7r = J" 00 dx (4>~ 2 {x) — A^ 2 ). In the asymptotic regime 



x < -£l (x > £r), we have X\(x) ~ A 



L(R) • 



and 



X n (x) ~ x/A L{R) + sgn(a;)7_AL(R), with 7_ = (tr - 
7l)/2. The solution to ((3|) can be written as 

X{x; k x ) = CMX^x; fc 2 ) + Cu(k x )X u {x; fc 2 ) (7) 

with X a (x; k 2 x ) ee J2n=o k 2 x n X ( a\x) for a = I, II and 

Xi°\x) =X a (x), (8) 



X^(x)=-4>(x) 



x dx 1 rx 



o <t> 2 (x')J 



d/^l")^- 1 '^"), 



(9) 
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for n > 1, where C a (k x ) are coefficients. Equating 
dU with (J7]) in the asymptotic regions £l,r "C \x\ <C 
\k x \~ 1 , we obtain lim fe ^ i(fc) = 2A- L A K /(A 2 ^ + A\) and 

lim fc ^oKfc) = ( A l^ A r)/( A 1 + A r)- We can then arrive 
at © with T{k) = \t(k)\ 2 and R(k) = \r(k)\ 2 . 

In the derivation of ([5]), the crucial facts were (i) that 
the WF of the excitation in the low-momentum limit is 
proportional to the condensate WF and (ii) that the heal- 
ing length for the condensate is finite so that the con- 
densate WF recovers the value in the spatially uniform 
system from the depleted value around the potential over 
a finite transient region. 

The Hamiltonian (JTJ is invariant under gauge trans- 
formation and spin space rotation, since V ex t is coupled 
only to the local density (i.e., symmetry-preserving po- 
tential). When continuous symmetry is spontaneously 
broken, the ground state has infinite degeneracy. Even in 
the spatially inhomogeneous system, the WF of the cor- 
responding excitation at E — > describes the NG mode. 
The fact that the spin wave is a NG mode is therefore 
linked to (i). The property (i) is unusual, since the WF 
obeying the Schrodinger equation equals the order pa- 
rameter of an interacting Bose gas. 

As for (ii), the finite healing length is characteristic 
of the interacting Bose gas; this interaction yields finite 
compressibility. When the compressibility is finite and 
(i) holds, the WF of the NG mode at the low-momentum 
limit as well as the condensate WF persist with finite 
and nonzero amplitudes near the potential barrier. This 
property of the WF leads to the nonzero transmission of 
the NG mode in the low-momentum limit. 

Scattering cross section. — We investigate the scat- 
tering cross section of the transverse spin wave from 
the spherical potential V eK t( r ) — V( r )> which satisfies 
lim r _ i-o r 2 V(r) = and linv-^oo V(r) = 0. In the ground 
state, the condensate WF is also spherically symmet- 
ric: </>(r) = 4>(r). The GP equation in dimensionless 
form is reduced to H r <t>(r) = with H r = -[(d 2 /dr 2 ) + 
(2/r)(d/dr)]/2 + V(r) and V(r) = V(r) + 4> 2 {r) - 1. 
<po (r) can be expanded with respect to the partial waves 
<^o( r ) = Ya^o Pi (cos 8) Ri(r; k) with the Legendre poly- 
nomial Pi (x) and the angle 8 between r and the incident 
momentum. The radial WF i?;(r;fc) satisfies 



rtr 



1(1 + 1) 



2r 2 



Ri(r-k) = —Ri(r;k), 



(10) 



with E = k 2 /2. It behaves as Ri(r;k) oc r l for r ~ 
and Ri(r; k) oc ji(kr) — tanSi(k)ni(kr) for r — > oo. Here, 
ji(kr) and ni(kr) respectively denote the spherical Bessel 
and Neumann functions, and 5i(k) denotes the phase shift 
of the partial wave. 

If V(r) were a generic potential, the phase shifts would 
behave as 8i(k) oc k 2l+1 for small k and the cross section 
cr(fc) = 47rfc~ 2 X^o(2' + 1) sm2 8i(k) would be given by 



for small k. However, in the present case, even when 
V(r) is generic, V(r) reflects the specific properties of 
the condensate WF. In the following few paragraphs, we 
show that this aspect has a significant influence on the 
s-wave scattering behavior, leading, at small k, to 



tan^o(fc) oc fc 3 , c(fc) oc fc 4 . 



(12) 



The derivation of (TT21 relies on the following two relations 
which correspond to (i) and (ii) noted previously: 



lim Ro(r; k) oc <f>(r), 

k— >0 

lim r m ((j)(r) - 1) = 0, 



(13) 
(14) 



for m > 0. In what follows, we fix the overall normaliza- 
tion factor of Ro(r; fc) such that 



R (r = 0;fc) =0(r = 0). 



(15) 



Let us expand Ro(r;k) as r^ 1 J2^ =0 fc 2n M^"^(r) with 
respect to fc 2 . It suffices to retain the first two terms 
tif°'(r) and u^(r) to prove (fT"2)) . u^(r) and vS l '(r) re- 
spectively satisfy 



H' r u^(r) = 0, H' r vS l >(r) = «W(r)/2, (16) 

where li' r = —\-^ps + V(r). A set of linearly indepen- 
dent solutions to the first equation of (TT^l) is given by 
ui(r) = r<p(r) and u u (r) = m(r) J^ (ui(r , ))~ 2 dr'. uj(r) 
is obtained by noting (TT3|) . The condition (fT"5|) yields 
«(°)(r) = Ui(r), since the second solution behaves as 
Un(r) ~ r° near r = 0. A special solution to the sec- 
ond equation of (|16[) is, on the other hand, given by 

vS^(r)=Ui(r) I ui(r')u\i(r')dr' —u\\(r) j (u\(r ')) 2 dr '. 
Jo Jo 

(17) 

With use of (|T7|) and u^°'(r) — Ui(r), we examine the 
asymptotic behavior of 

R (r;k) ^u^\r)/r + k 2 u { - 1 \r)/r + 0(k 4 ), (18) 

at r — > oo. (fTl)) yields u^°'(r)/r = u\(r)/r ~ 1 and 
u w (r)/r ~ -r 2 /6 + B + C/r at r -> oo. Here, B and C 
are given by 

/•OO />oo 

B= [ui(r')uii(r')+r']dr', C= [(wi(r')) 2 -r' 2 ]dr' '. 
Jo Jo 

(19) 

The asymptotic form of Ro(r;k) then reads 



> 0I Wm - „(o)/ 



Ro(r;k) = (l + Bk 2 ) 1 



fc 2 r 2 
~6~ 



Cfc 3 



-1 
kr 



cr(fc) ~ fc 2 sin 2 S (k) ~ const. (> 0) 



(11) 



(20) 

Since jo(kr) and n^(kr) are respectively given by 
jo(kr) ~ 1 — fc 2 r 2 /6 and no(fcr) ~ — l/(fcr) for kr -C 1, 
the phase shift of the s-wave (I = 0) can be written as 
tan<5 — Cfc 3 /(1 + Bk 2 ) oc fc 3 . flTJ) has been proved. 



4 



There are known examples in which a(k) oc k 4 holds for 
the wave with E oc k; these examples include sound, light 
(Rayleigh scattering) , and the Bogoliubov mode [5j, |6| of 
a BEC. Our result a(k) oc k 4 for low-energy scattering of 
the spin mode with E oc k 2 is rather unfamiliar compared 
to (jlll) and is regarded as an anomalous power law. 

Discussion. — The transmission coefficient T of the Bo- 
goliubov mode in the spin-1 BEC is the same as that in 
the scalar BEC; we here recall this coefficient of the Bo- 
goliubov mode in the low-energy limit 0, [2l| given by 



4 A h A R cos 6>l cos R 
hm 1 = — — — , 

fc->o (Al cos 6l + j4r cos Wr y 



(21) 



where 0r is the refracted angle and others are defined 
previously. Comparing ([5]) and (121 1) with the transmis- 
sion coefficient of classical waves, we see that the mod- 
ulus of magnetization M„ and that of the condensate 
WF A v both play the role of "impedance" . We can then 
extract the following common feature of both low-lying 
modes: (a) In a system whose continuous symmetry is 
spontaneously broken, impedance of the corresponding 
NG mode in the transmission coefficient is given by the 
amplitude of its order parameter (e.g., |(^)| or |M|). We 
omitted the subscript j = +1 of As for the low- 
energy scattering from a spherical potential, we previ- 
ously noted another unexpected feature common to the 
two low-lying excitations: (b) the s-wave scattering cross 
section is proportional to the fourth power of the momen- 
tum. Table U summarizes the scattering properties of the 
two NG modes. The derivations of (a) and (b) both 
rely on (i) and (ii), and are straightforward and general 
without use of any other property specific to the system. 
We thus anticipate that (a) and (b) both hold as general 
properties of NG modes in the presence of symmetry- 
preserving potentials 22 1. 



is negligible, the following two s-wave scattering cross 
sections will be observed: the first is the Rayleigh type 
(fl~2f for a doped nonmagnetic impurity (e.g., a closed 
shell atom) or for a defect like a mechanical depression of 
the magnetization. The second is a usual type (fTTj) for a 
doped magnetic impurity (e.g., an atom classified as tran- 
sition metals or rare earth elements). Propagation of the 
spin waves can be manipulated with impurities and de- 
fects adequately disposed. It was recently reported that 
electric signals in a magnetic insulator are transported 
by converting the electric current to a spin wave [23j j . 
If (fT2|) holds and there are no magnetic impurities, the 
scattering of spin waves from nonmagnetic impurities and 
inhomogeneities is suppressed for lower momentum, and 
electric signals could be transported over longer distances 
at lower energy. 

Note that, in ultracold gases, the diffraction and the 
scattering problem of excitations from a spherical poten- 
tial are now possible to investigate if one uses a single 



24] 



ion trapped inside a BEC 

Conclusion. — We constructed the simplest model of 
an anomalous tunneling phenomenon for an excitation 
through a generic localized potential, using a transverse 
spin wave in the ferromagnetic phase of a spin-1 BEC. 
Interestingly, studies of anomalous tunneling will connect 
to many systems whose continuous symmetry is sponta- 
neously broken. Our present study will be a minimal 
model for these studies. 
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TABLE I: Scattering properties of NG modes in the ferro- 
magnetic phase of the spin-1 spinor BEC 



type of scattering 

NG mode NG mode refraction impedance cross section 

Bogoliubov I Snell's law 1(^)1 o"o oc k 4 

transverse spin II no |M| oq oak 4 



It was demonstrated that spin waves in magnetic sys- 
tems transmit a barrier induced by the inhomogeneous 
magnetic field [ItJ ■ This high tunneling probability is due 
to the nonlocal magnetic dipole interaction of spins. We 
note that this is different from our situations; we propose 
a completely different mechanism of high tunneling prob- 
ability through the nonmagnetic barrier. It stems from 
the degeneracy of the ground state induced by the spon- 
taneous breaking of the continuous symmetry. In mag- 
netic systems whose nonlocal magnetic dipole interaction 
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